Frontiers in Molecular Biosciences

Background: Citrin deficiency, caused by biallelic pathogenic variants in SLC25A13, must be identified early to prevent serious complications such as hyperammonemia and liver failure. However, clinical diagnosis is often delayed due to its nonspecific presentation and limited sensitivity of amino acid-based newborn screening methods. Although genome-based evaluations are being investigated to address these issues, concerns about their cost, turnaround time, variant interpretation ability, and data handling highlight the need for a more practical yet reliable alternative. We investigated the feasibility of applying proteomic approach on dried blood spots (DBS), which are routinely used in newborn screening. Methods: We performed untargeted liquid chromatography-tandem mass spectrometry to analyze the proteome of DBS using a previously developed "non-targeted analysis of non-specifically DBS-absorbed proteins" (NANDA) workflow. SLC25A13 protein abundance was quantified in individuals with biallelic loss-of-function mutations, compound loss-of-function/missense mutations, and heterozygous carriers; this was also evaluated in healthy and diseased controls representing relevant differential diagnoses. To leverage proteomic information, we derived a multivariate proteomic signature using feature selection and evaluated its performance with leave-one-out cross-validation. Biological relevance was assessed by enrichment analysis, and complementary transcriptomics was performed using RNA sequencing. Results: A total of 7,474 proteins, including SLC25A13, were consistently detected in DBS. SLC25A13 was undetectable in individuals with biallelic loss-of-function mutations. However, individuals with compound loss-of-function/missense genotypes showed reduced but measurable SLC25A13 levels, comparable to those observed in heterozygous carriers. In contrast, a compact 15-protein signature accurately identified individuals with compound loss-of-function/missense genotypes (AUC, 0.99; sensitivity, 1.00; specificity, 0.95). The signature was enriched for Ca2+-response, and transcriptomics showed downregulation of genes related to multimodal ion channels in affected individuals compared to controls. Conclusions: DBS-based proteomic profiling may assist in the diagnosis of citrin deficiency through SLC25A13-quantification and a biologically plausible multivariate signature. More broadly, this strategy offers a promising new diagnostic layer for protein disorders, providing a proteomic readout in a clinically practical DBS format with potential utility for future diagnostic and screening applications.

Tyrosine Kinase inhibitors (TKIs) are widely used as effective chemotherapeutic agents for treating patients with EGFR-mutated NSCLC. Unfortunately, after treatment, patients eventually develop resistance to TKI therapy. The most common resistance mechanism for the TKI Osimertinib is the overexpression of the MET Proto-Oncogene, Receptor Tyrosine Kinase (MET). We previously demonstrated that miR-411-5p A-to-I edited at position 5 (miR-411ed) can directly target MET in A549 and H1299 cells. MiR-411ed in combination with Osimertinib reduced cell proliferation in two TKI resistant EGFR-mutated cell lines: HCC827R and PC9R. MiR-411ed did not downregulate MET expression in HCC827R, suggesting an alternative mechanism for TKI response. In this study, we aim to identify the mechanism of miR-411ed TKI response using a multi-omics approach of RNAseq and protein mass spectrometry. In our cellular model, we identified miR-411ed affected genes independent of MET activity, resulting in 211 genes (RNAseq) and 36 proteins (proteomics). Pathway analysis identified an increase in interferon signaling for RNAseq and combined omics, and a decrease in ERK/MAPK signaling in proteomics. Using the IsoTar target prediction tool, we identified STAT3 as a key regulator and confirmed STAT3 protein downregulation upon transfection with miR-411ed. We further investigated the effect of miR-411ed in vivo, observing a reduction in tumor size with miR-411ed in combination with Osimertinib but not with miR-411ed or Osimertinib treatment alone, confirming the effectiveness of miR-411ed in TKI response.

Proteins with intrinsically disordered regions (IDRs) migrate at a higher apparent molecular weight in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) complicating their analysis and identification. Here, we investigate the sequence determinants of the hypomobility of IDRs using a series of synthetic low complexity domains. We find that negative charge increases the apparent molecular weight, but neutral polar tracts also have abnormally slow migration. Positive charge and hydrophobic residues decrease the apparent molecular weight, although lysine residues show a biphasic effect with decreased migration at high fractional contents. Combinations of residues show that different sequence contributions to the apparent molecular weight are not additive. The results can be rationalized by the protein-decorated micelle model by considering both SDS binding and the compaction of protein SDS-complexes.

Breast cancer (BC) remains a leading cause of morbidity and mortality worldwide. Early detection remains the most effective strategy for improving prognosis. We explored the urinary extracellular vesicle (uEV) proteome for changes linked to BC which could also be potential biomarkers. Urine samples were collected from 20 participants across four groups (n = 5 each): newly diagnosed BC patients, benign breast disease (BBD) patients, individuals with breast cancer symptoms (symptom control, SC), and age-matched healthy controls (HC). EVs were isolated using size exclusion chromatography and extracted proteins were analysed using a GeLC proteomic approach. Proteins were identified and quantified using Proteome Discoverer and further analysed using MetaboAnalystR, Funrich and Metascape. A total of 256 proteins were identified from the uEV preparations. BC comparisons with BBD, SC and HC identified 7 proteins differentially expressed proteins (DEP); SERPINB1 -- Serpin family B member 1, LCN1 -- Lipocalin 1, SIRPA -- Signal regulatory protein alpha, ACTB -- Actin, beta, YWHAZ --Tryptophan 5-monooxygenase activation protein zeta, Ig JCHAIN and APOA1 -- Apolipoprotein A1. Receiver Operator Characteristic (ROC) curve assessments suggested that each DEP protein had an area under the curve (AUC) of > 0.8. These findings highlight EV-derived proteins as promising non-invasive biomarkers for breast cancer detection, warranting further validation in larger cohorts.

Staphylococcus aureus (S. aureus) is an opportunistic pathogen that is a global health concern for its ability to cause a wide spectrum of clinical infections. Due to the emergence of resistance to commonly used antibiotics, there has been interest in exploring the use of antimicrobial peptides to treat S. aureus infections. However, changes in the lipid composition of the lipid bilayer membrane can alter the activity of peptides, and S. aureus is able to induce variations in lipid composition in response to environmental stress. Here, we explore how the main lipid components in S. aureus are altered when exposed to LL-37, a human cathelicidin involved in primary immune response, and ATRA-1, a short antimicrobial peptide derived from the snake Naja atra venom. A lipidomic study is conducted through HPLC-MS-MS (LC-ESI-MS/MS) to quantify phosphatidylglycerol, cardiolipin, lysyl-phosphatidylglycerol, monogalacto- and digalacto-diacylglycerol, and carotenoids. In addition, menaquinones, responsible for electron transport during oxidative phosphorylation, were also quantified. Biophysical properties such as membrane electric surface potential and lipid packing were assessed. We find that lipid adaptation is specific to the type of antimicrobial peptide, where ATRA-1 mainly induces changes in the electric surface potential through variations in Lysyl-PG, while exposure to LL-37 changes carotenoid levels, inducing an increase in membrane rigidity as measured by FTIR. In addition, both peptides induce a reduction in menaquinone and DGDG levels. These findings highlight the role of membrane lipid remodeling as a peptide-specific response mechanism in S. aureus, with implications for the development of AMP-based therapies. HighlightsO_LIStaphylococcus aureus responds through shifts in lipid composition and membrane biophysical properties to exposure to the antimicrobial peptides LL-37 and ATRA-1. C_LIO_LIBoth LL-37 and ATRA-1 lead to shifts in the glycolipids MGDG and DGDG; two lipids involved in regulating negative membrane curvature stress and responsible for shifting resistance to antimicrobial peptide activity in Staphylococcus aureus. C_LIO_LILL-37 treatment leads to an overall reduction in carotenoid content in Staphylococcus aureus, including the carotenoid end-product staphyloxanthin and the precursor 4,4-diaponeurosporenoic acid. Both lipids regulate membrane biophysical properties and protect Staphylococcus aureus from oxidative stress. C_LIO_LIBoth LL-37 and ATRA-1 lead to a reduction in menaquinone levels, which are involved in the electron transport chain during oxidative phosphorylation. Reduction in these menaquinones have been associated to the formation of small colony variants that are often observed in chronic Staphylococcus aureus infections. C_LI

Staphylococcus aureus (S. aureus) is a clinically relevant pathogen capable of adapting its membrane composition in response to environmental stress. In this adaptive process, bacterial carotenoids play a crucial role. Although staphyloxanthin (STX) is the main carotenoid produced by the bacterium, S. aureus also synthesizes other pigmented intermediates that play an unknown role in regulating membrane biophysical properties. In this study, we purified 4,4-diaponeurosporenoic acid (4,4'-DNPA) from S. aureus carotenoid extracts and evaluated its effect on the thermotropic and biophysical properties of representative membrane models. The highly rigid triterpenoid 4,4'-DNPA is one of the last precursors in the biosynthesis of STX and is found in high concentrations in the stationary phase of S. aureus. Phase transition temperatures were determined using infrared spectroscopy, while interfacial hydration and hydrophobic core dynamics were investigated using fluorescence spectroscopy through Laurdan generalized polarization and DPH anisotropy. The results show that 4,4'-DNPA increases the main phase transition temperature of lipid bilayers in a concentration-dependent manner. This is in contrast to STX that decreases the transition temperature. This difference is consistent with the additional fatty acid present in STX that changes its effect on the phase behavior. Furthermore, 4,4'-DNPA reduced the interfacial hydration levels and restricted hydrophobic-core dynamics at higher concentrations, consistent with increased molecular order and stability. 4,4'-DNPA therefore complements STX in increasing membrane order and lipid packing. These findings support the notion that the production of bacterial carotenoids functions as a biophysical regulatory mechanism of lipid packing in S. aureus membranes.

The HIV-1 Vpu protein aids viral adaptation by influencing host cell pathways via protein interactions. While Vpu is mainly found in plasma and endomembranes, we recently discovered a soluble form that forms a stable, equimolar complex with Ca2+-bound calmodulin (Ca2+-CaM), potentially affecting Vpus cellular trafficking. Here, to determine the binding affinity and identify regions of soluble Vpu involved in CaM binding, we used ensemble Forster Resonance Energy Transfer (eFRET). We tested Cy3-labeled full-length (FL) Vpu, a C-terminal fragment (helices 2 and 3), and a Cy3-labeled FL Vpu V22A/W23Y mutant with substitutions in Vpus helix 1. All Vpus variants were labeled at residue L42C, and Ca2+-CaM was tagged with Cy5 at residue S39C. eFRET analysis of 100 nM Cy3-Vpu variants mixed with Cy5-Ca2+-CaM (in the range 100-600 nM) revealed dissociation constants (Kd) and binding energies ({Delta}G) for heterocomplexes. FL Vpu-Ca2+-CaM showed high stability (Kd [~]40 nM,{Delta} G [~]10.1 kcal/mol), while the truncated C-terminal region and V22A/W23Y mutant formed less stable complexes with Ca2+-CaM (Kd[~]200 nM and 800 nM,{Delta} G [~]9 kcal/mol and [~]8.3 kcal/mol). This, a binding hot spot in Vpus CaM-binding motif in helix 1 was identified, which may control the stability of Vpu-Ca2+-CaM complex and Vpus insertion in the membrane: We hypothesize that upon delivery to the membrane, the hydrophobic helix 1 of Vpu dissociates from Ca2+-CaM and inserts in the lipid bilayer; thereafter, CaM dissociates from Vpu facilitated by the reduced Vpu-Ca2+-CaM complex stability. The findings from this study advance our understanding of HIV-1 Vpu interactions with cellular components and may aid the development of antivirals.

Inflammasomes are cytosolic innate immune sensors that, once activated by a pathogenic threat, lead to activation of the inflammatory Caspase-1. Inflammasome activation and its consequences have been studied extensively in myeloid cells and in overexpression systems. Recent studies have identified cell type specific effects that are not fully explained by the known cleavage targets of Caspase-1. Here, we identified targets of caspase cleavage using mass spectrometry in primary intestinal epithelial cells by specifically activating the NAIP-NLRC4 inflammasome. We have taken an unbiased approach and developed a novel method for analyzing mass spectrometry data for evidence of caspase activity. Our approach can also be applied to existing proteomic datasets to establish the presence of caspase activity under various biological conditions. These results lay the groundwork for future studies on mechanisms of caspase-induced processes such as intestinal epithelial cell extrusion.

Although several ancillary tests are available in limited laboratories, diagnosis of microphthalmia (MiT)/TFE family translocation renal cell carcinoma (tRCC) could be challenging due to diverse and overlapping tumor morphology and the lack of reliable biomarkers. GPNMB has been recently identified as a diagnostic marker for various renal neoplasms with FLCN/TSC/mTOR-TFE alterations. However, the sensitivity and specificity of GPNMB immunostain are suboptimal and the result interpretation in ambiguous cases could be difficult. To search additional biomarkers that could improve the screening sensitivity and predict genetic aberrations in FLCN/TSC/mTOR-TFE pathway in renal tumors, we performed bioinformatic analysis of publicly available cancer databases and found GPR143, a transmembrane protein regulated by MiT transcription factors, was highly expressed in a subset of renal cell carcinomas (RCCs). In two the Cancer Genome Atlas (TCGA) kidney cancer cohorts, RCCs with high levels of GPR143 expression were enriched for renal neoplasms with FLCN/TSC/mTOR-TFE alterations. Similar to GPNMB labeling, GPR143 immunostain was positive in the majority of tRCC cases and renal tumors with FLCN/TSC/mTOR alterations, suggesting that GPR143 could function as another surrogate marker for FLCN/TSC/mTOR-TFE alterations in certain renal tumors. Interestingly, despite the concordant GPR143 and GPNMB immunoreactivity in most renal neoplasms with FLCN/TSC/mTOR-TFE alterations, diffuse GPR143 immunostain was observed in some cases with negative or focal GPNMB labeling. Taken together, our results indicate GPR143 could serve as a useful adjunct marker to improve the sensitivity for screening renal tumors with FLCN/TSC/mTOR-TFE alterations.

Lipoprotein composition is altered in sepsis, and supplementation with high-density lipoproteins has been reported to improve outcomes in experimental settings. In this study, we aimed to investigate the nature and inter-individual variability in the lipoprotein proteome to inform risk stratification and opportunities for precision medicine approaches. In a large proteomic dataset including 1134 patients (1781 samples) with sepsis and 149 healthy volunteers, we analysed 18 protein components of lipoproteins. We characterise heterogeneity of the lipoprotein proteome, defining three step-wise sub-phenotypes associated with increasing disease severity, one close to health, then an early phase patient group showing increased abundance of proteins that integrate HDL under inflammatory conditions (SAA1 and SAA2), then a group with decreased abundance of proteins that are components of HDL under healthy conditions that was associated with higher organ failure intensity (SOFA score) and increased mortality. We developed and externally validated a quantitative score reflective of lipoproteins alterations in sepsis, and machine learning predictive models to predict the LP class, advancing future individualised lipoproteins-based therapeutics in sepsis.

Intercalation of small molecules between DNA base pairs affects DNA conformation, disrupting essential cellular processes including replication, transcription, and repair. We investigated conformational changes in 18-mer DNA upon intercalation of doxorubicin, SYBR Gold and YOYO-1 using extensive MD simulations. Two main patterns for the intercalation were identified: RISE-type intercalation occurs between adjacent base pairs and extends the DNA helix with decreased twist angles, while OPEN-type intercalation proceeds through base-pair opening without significant DNA extension. Kinetic analysis revealed that association rates for intercalation followed the order: first YO-moiety (mono-intercalation) > SYBR Gold > doxorubicin > YOYO-1 (bis-intercalation). Free energy landscape showed that forces at DNA termini reached up to 117 pN during stretching. Notably, base pairs adjacent to intercalators were protected from strand separation, accompanied by additional helical unwinding. MM-PBSA/GBSA analysis revealed that the driving force for intercalation is the stacking energy, and the binding affinity was highest for minor groove binding. Persistence length decreased with single molecule binding but recovered with two molecules due to their electrostatic repulsion. Mechanical properties of intercalated DNA showed position-dependence, demonstrating that multiple intercalation modes coexist in solution. The heterogeneous nature of intercalation explains why experimental measurements reflect ensemble averages rather than single binding configurations.

This study investigates the interaction between the cationic antimicrobial peptide protamine and bacterial porin OmpF (E. coli) at the single-molecule level. Using high-resolution conductance measurements in planar lipid bilayers, strong voltage- and concentration-dependent ion current blockages with OmpF, indicating significant protamine binding were observed. Further analysis revealed that peptide length influences binding kinetics, with longer peptides showing reduced affinity and slower exchange rates. These findings demonstrate that OmpF is a tractable model for studying cationic peptide-channel interactions and translocation mechanisms relevant to antimicrobial action.

Biomolecular condensates that form via liquid-liquid phase separation (LLPS) of, most prominently, intrinsically disordered proteins (IDPs) are ubiquitous in eukaryotic cells and responsible for regulating a plethora of biological functions. Amongst these, they contribute to regulating cell motility, either individually within an extracellular matrix or collectively within confluent epithelial tissue. In this computational study we focus on the latter with the aim of investigating whether the mutual exertion of mechanical forces during collective migration in an epithelium can principally trigger cytoplasmatic LLPS. Since present models for confluent epithelial motility have so far only considered cells that are devoid of phase separating (protein) solutes, we extend a common multiphase approach for 2D cell motility with a mixing contribution including any number of protein solutes. Our model considers the phase behavior in both intracellular and extracellular regions and determines to what extend the membrane is permeated by the solutes under the influence of mechanical and osmotic forces. Our initial calculations unlock a very rich behavior involving formation and dissolution of condensates during migration, as well as an impact of LLPS on the very nature of the motility itself, through feedback mechanisms which may bear biological relevance.

Coarse-grained structure-based models (CG-SBMs; or G[o] models) are simplified potential energy functions of biomolecules or biomolecular complexes that encode their structure. Molecular dynamics simulations of such SBMs have been successfully used to study long time-scale dynamics such as protein and RNA folding, and large conformational transitions of biomolecular complexes. SBMs have several advantages: (1) Their MD simulations are computationally inexpensive, making extensive sampling easily accessible to many researchers. (2) They are easy to modify and can be adapted for the specific biomolecular problem that needs to be investigated. However, the force-fields of SBMs are not usually included in commonly used biomolecular simulation packages resulting in a barrier to their use. Here, we present SuBMIT (Structure Based Models Input Toolkit; https://github.com/sglabncbs/submit), a toolkit for generating coarse-grained SBM input files for performing MD simulations with GROMACS and OpenMM/OpenSMOG. Simulations whose input files can be generated using the different flavors of CG-SBMs present in SuBMIT include the folding and conformational ensembles of proteins with intrinsically disordered regions, 3D-domain-swapping in proteins and the dynamics of RNA-protein assemblies (e.g., simple RNA viruses).

Surface factors that contribute to the virulence of Staphylococcus aureus have become therapeutic targets in the treatment of illness associated with this bacterium. Staphylococcal protein A (SpA) is a well-known contributor to S. aureus toxicity and virulence, although relatively little is known about protein A and how its biological function has evolved. SpA is displayed on the surface of the bacterium and contains 5 nearly identical helical ({approx} 60 aa) domains that bind antibodies with high affinity (Kd {approx} 10 nM). The folding free energies of only domains E and B have been determined. In this study we used intrinsic fluorescence detected denaturation to measure the folding thermodynamics of each domain in isolation and in the native multidomain context using a construct that includes the N-terminal half of the mature protein (SpA-N). We also constructed a series of proteins with 1 to 5 repeats of B domain, linked exactly as the five domains of WT SpA are linked. We used nearest neighbor thermodynamic models to explicitly demonstrate that the domains in B domain repeat proteins fold independently. We also showed that the domains in SpA-N fold independently by comparing the folding free energies of domains in isolation and in their multidomain context. Previous dynamic NMR experiments detected highly flexible linkers between domains in 5B, suggesting that the domains of SpA are structurally independent, which is likely responsible for the lack of thermodynamic coupling. Our results also showed a steep increase in domain stability from the N-to C-terminus in SpA-N, from 0.97 {+/-} 0.05 to 5.57 {+/-} 0.28 kcal/mol. We hypothesize that this stability gradient is related to efficient secretion of protein A.

The importance of human aldehyde oxidase (hAOX1) has increased over the last decades due to its involvement in drug metabolism. Inhibition studies concerning hAOX1 are extensive and a common reducing agent, dithiothreitol (DTT), was recently found to inactivate the enzyme. However, in previous crystallographic studies of hAOX1, DTT was found to be essential for crystallization. To surpass this concern another reducing agent used in crystallization trials. Using tris(2-carboxyethyl)phosphine (TCEP), a sulphur-free reducing agent, it was possible to obtain well-ordered crystals from hAOX1 wild type and variant, hAOX1_6A, which diffracted beyond 2.3 [A]. Instead of the typical star-shaped crystals of hAOX1, at pH 4.7, plates are obtained in the orthorhombic space group (P22121) with two molecules in the asymmetric unit. Activity assays with the enzyme incubated with both reducing agents show that contrary to DTT, TCEP does not lead to irreversible inactivation of the enzyme. The replacement of DTT with TCEP in crystallization of hAOX1 provides a strategy to circumvent enzyme inactivation during crystallographic studies, allowing future applications of new assays, such as time-resolved crystallography.

1Long QT syndrome Type 2 (LQT2) is a genetic disorder caused by missense mutations in the KCNH2 gene that encodes the potassium channel KV11.1. Previous studies have shown that most KV11.1 missense mutations with loss-of-function phenotypes result from impaired trafficking from the endoplasmic reticulum to the plasma membrane. To investigate the molecular basis of these defects, we used molecular dynamics simulations to analyze two sets of disease-associated missense mutations: those that suppress and those that maintain normal channel trafficking. We focused initially on the conformational and dynamics differences between wild-type and several mutants of KV11.1 via molecular dynamics simulations when two K+ were placed in the selectivity filter (SF). Our study reveals that missense mutations in the S4 helix allosterically disrupt the selectivity filter, a critical determinant for proper channel trafficking. Trafficking-competent variants largely retained a wild-type selectivity filter structure, whereas trafficking-deficient mutants exhibited pronounced structural perturbations in this region. These findings suggest that certain LQT2-associated missense mutations in KCNH2 impair channel trafficking by compromising the structural integrity of the selectivity filter. We additionally found that second-site variants Y652C in the drug binding vestibule can correct structural defects associated with some mistrafficking variants.

The Na+/K+-ATPase (NKA) regulates ion balance in the kidney and influences cellular processes like proliferation and apoptosis through its signal transduction. The endogenous ligand 20-Hydroxyeicosatetraenoic acid (20-HETE) contributes to inflammation and fibrosis in chronic kidney disease (CKD) and inhibits NKA activity in renal tubules. However, the molecular mechanism of this interaction remains unclear. In this study, we used in-silico approach to investigate the potential interaction between 20-HETE and NKA. Various ligands, including known NKA ligands such as cardiotonic steroids (CTS), 20-HETE, and negative controls, were docked using rigid and Induced Fit Docking to predict the affinity of the ligands toward NKA. Binding free energy calculations with the Prime Molecular mechanics with generalized Born and surface area (Prime MM/GBSA) tools were used to confirm the involvement of key amino acids in ligand-receptor interactions. The docking analyses revealed that 20-HETE exhibited a binding affinity comparable to negative control, with some differences between rigid and induced fit docking. Binding free energy data highlighted key amino acids in the 20-HETE and NKA interaction. Interaction fingerprint and mutations such as Ala330Gly and Val329Ala significantly reduced binding free energy, while Thr804Ala showed a notable decrease, underscoring the potential importance of these amino acids in ligand stabilization. These findings provide computational evidence supporting potential direct interaction between 20-HETE and NKA and identify candidate residues for future experimental validation.

G-quadruplexes (G4s) are critical nucleic acid secondary structures that play pivotal roles in regulating gene expression. In this study, we conducted a proteome-wide in silico analysis across multiple viruses causing hemorrhagic fevers to identify candidate proteins containing a conserved G4-binding motif. Four peptides belonging to Marburg, Ebola, Hantaan and Yellow fever viruses were shown to bind to G4 in vitro. We selected the NS3 protease domain of Yellow Fever virus for further validation. Biochemical assays demonstrated that the NS3 protease domain binds G4 structures with high specificity and affinity, particularly favoring the parallel conformation. Molecular docking and simulations further revealed that the NS3 protease domain interacts with the terminal G-tetrads and loop regions of G4 via key residues, including PHE40, adopting an insertion and stacking composite binding mode. These findings expand our understanding of virus - G4 interactions and offer novel potential targets for G4-based antiviral strategies. Bullet points- We screened viruses causing hemorrhagic fevers for potential G4-binding peptides. - Four peptides belonging to Marburg, Ebola, Hantaan and Yellow fever viruses were shown to bind to G4 in vitro. - Biochemical assays demonstrated that the NS3 protease domain of YFV binds G4 structures with high specificity and affinity.

In any field, unquestioningly accepting artificial intelligence (AI) results should be considered bad practise. Here, we devised a comparative modelling-based strategy for validating protein structures that exploits the well-known observation that protein folds are far more conserved than protein sequences. We identify proteins with a similar fold to the AlphaFold-generated query protein and determine their structural alignment to the query. The hypothesis is that if the sequence alignment coincides with the structural alignment, then the structure is validated. The strategy is implemented on a helix-by-helix and strand-by-strand basis using a multi-template pairwise local profile alignment method that works well into the twilight zone. The method is illustrated by application to the transmembrane transporter PEPT1, for which the structure is known, and the S-deacylases ABHD13 and ABHD16A, for which only AI-generated models exist. ABHD16A is particularly challenging because a sequence alignment search with BLASTp does not reveal any structural homologues and therefore requires work with extremely remote homologues; however, both models are validated through this strategy and are stable during classical molecular dynamics simulations. The ability of the strategy to identify errors is assessed with reference to misaligned ABHD13 models and misfolded decoy proteins.